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Abstract 

We introduce a way of implementing Wilson renormalization within the con- 
text of the theory of effective Hamiltonians. Our renormalization scheme 
involves manipulations at the level of the generalized G-matrix and is inde- 
pendent of any specific kinematics. We show how to calculate the beta func- 
tion within this context and exhibit our method using simple scale-invariant 
quantum mechanical systems. 
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The theory of effective Hamiltonians and operators in many-body physics has a long 
and elaborate history This framework primarily involves manipulation of operators 
and is independent of the chosen kinematics. Our present effort is aimed at incorporating 
the concepts of renormalization within the context of the theory of effective Hamiltonians. 
Recently, other approaches for Hamiltonian renormalization have appeared both in field 
theory and quantum mechanics ||.[| . 

Our general philosophy towards renormalization is inspired by the original work of Wilson 
f|. Wilson's approach involved integrating out degrees of freedom above a certain momen- 
tum range to arrive at an effective action. This concept has primarily been implemented in 
the path integral formulation of a given problem. In recent years there has been a renewed 
effort to implement renormalization within a Hamiltonian formulation 0-0], as this may 
lead to a better understanding of key issues in strong interaction physics. 

In this Letter we present a new non-perturbative scheme for renormalization which 
utilizes some of the more recent developments in the theory of effective Hamiltonians for 
many-body systems. Although our approach is quite general, we illustrate the methods with 
simple one-body quantum mechanical problems. In this way, we can exhibit the efficacy of 
our ideas without being overwhelmed by technical complications. 

Consider the eigenvalue problem 

H\* i ) = E i \%) z = l,2,...,n, (I) 

with n possibly infinite. There are very few Hamiltonians for which Eq. ([!]) is exactly 
solvable. In addition, the usual perturbation theory may turn out to be inadequate, and 
a large number of degrees of freedom associated with a many-body system may prevent a 
straightforward analysis of the problem. In realistic situations one would, in general, seek a 
subset of all the solutions of the above Hamiltonian. This naturally leads to the concept of 
an effective Hamiltonian. 

We shall now briefly describe our construction of the effective Hamiltonian. We can 
arbitrarily split up any Hamiltonian H as 
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H = H + V, 



(2) 



where Hq is exactly solvable. One motivation for doing this step is to provide a convenient 
basis to work with, defined by the set of eigenvectors of H : 



tf |/i>=£». 



(3) 



All matrices will be written with respect to this basis, unless otherwise indicated. 

We choose a model space M. which contains d basis vectors of H . The operators P and 
Q which project into and out of our model space, respectively, are given by 



p= E Im)H 

Q= E l^| = i -p. 

It is helpful to think of these projection operators as partitioned matrices : 



(4) 
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Consider, as in 0, a transformation of our Hamiltonian 

H = e~ s He s , 
|* i )=e- s |* i ) J 



(5) 



(6) 



where S is an operator to be determined shortly. It follows directly from this transformation 
that the eigenvalues of H are the same as the eigenvalues of the original Hamiltonian, i.e., 



H\%) = Ei\%). 

We will use the freedom in the choice of S to require that 

H{P\%)) = z = l,2,...,d. 
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(7) 



(8) 



As expressed in conventional applications to many-body problems, the goal of the effective 
Hamiltonian formalism is to construct an operator which acts only in a model space, yet 
gives us a subset of the exact eigenvalues of the full Hamiltonian. For our purposes here, 
we restate the goal as that of obtaining a subset of well-defined solutions of the eigenvalue 
problem for H. 

For this purpose, we choose the effective Hamiltonian as 



H eS = PHP, 



(9) 



which clearly acts only on the states in the model space (which may, in itself, be infinite 
dimensional). From Eqs. (S) and (RJ) it also follows that 



H f a(P\%))=E i (P\* i )). 



(10) 



The expression in Eq. (|9|) is therefore a consistent choice for H c g. 

The problem of finding H e g now reduces to one of finding an appropriate S. Following 
H|7|], we choose to obtain S such that S = QSP. This implies that in our chosen basis S 
must have the form 



S 



(11) 



where s is a (presently) arbitrary n — d by d matrix. 

By choosing 5* in this particular way, S n is zero for all n > 1. This leads to 



l + S 



with an analogous result for e s . 

If we started out with an arbitrary hermitian Hamiltonian 



fl 




V 





(12) 



H 



a 


b 




f 



\ 



(13) 



with a = and f — f\ then 



H = e~ s He s 
( 



a + bs 


b 


-s(a + bs) + fct + fs 


f-sb 



(14) 
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The effective Hamiltonian H e g therefore takes the form 

H eS = PHP = a + bs. (15) 
We shall now exhibit an iterative method to obtain H e g. Following 0] we define 



Z = H cS — u = a + bs — u>, 



(16) 



which is equal to the effective Hamiltonian up to the arbitrary additive constant u. Next 



Hn 
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(17) 



where Ap And Aq are the (diagonal) matrices containing the eigenvalues of Hq. Finally we 
introduce a generalized G-matrix defined as [6-7] 



G(u) = PVP + PVQ 
= PVP + PVQ- 



uo-QHQ 
1 



QVP 



:QVP + 



uj - QH Q 

+PVQ kwn QVQ 7TW7^ QVP + 

which, for the above conventions, can be written as 

G(uj) = (a-X P ) + b^—bl 
uj- f 



(18) 



(19) 



Eqs. (|H]) and ( |T%D can be solved iteratively to give Z and G(uj). From now on, we will 
denote G(uj) as G. One such iteration scheme is M 
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Z x = PHqP + G-luP, 

Zfl -i f~i j^i ry (~t ry ry 1^*7*7 V ^1 5 (^^) 

1 — trl — Or2^n-l — <-*3^n-2^n-l — • • • — Kj n -\Zj2^'i ' ' ' 



where 



H e s can finally be constructed from the above solution for Z. 

We note in passing that the generalized G-matrix may provide a leading approxima- 
tion to H c g. Within that approximation, our Eq. ([18|) bears resemblance to the effective 
Hamiltonian introduced in [[|. 

We shall now introduce the concept of renormalization within the above framework. We 
have seen above that the knowledge of the matrix G allows us to obtain Z, which is identical 
to if e ff up to an additive constant. In what follows, we shall therefore restrict our attention 
only to G. For the sake of convenience we choose to work in the momentum representation 
where the kinetic energy term in the Hamiltonian is diagonal. To introduce the concept 
of renormalization we shall focus our attention on the one-particle system. The formal 
generalization to a many-particle system would be straightforward. The matrix elements of 
G are here given by 

G kk , = (k\PVP\k') + 

+ [ dpdp'(k\PVQ\p)(p\ L^\p')( p '\QVP\k') + --- (22) 

Let us suppose that the potential V depends on a single coupling constant /x , which we 
shall call the bare coupling constant. It is clear from Eq. ( p2j) that the matrix element Gkk' 
will be a function of /io- The expression in Eq. (p2|) may, in general, require regularization 
due to the divergence arising from the integral. The regularization that we choose consists 
of introducing an ultraviolet cutoff A. The matrix element in Eq. ( p2| ) is now a function of 
the coupling constant /x and the cutoff A. At the end of the calculation we must remove 
the cutoff, i.e. we must take A to oo, which, as discussed above, may in general lead to 
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divergence. One way to avoid the divergence is to replace the coupling constant [j,q with a 
function of A, which we denote as /i(A), and then require that matrix element in Eq. ( |22|) 
remain finite and independent of the cutoff as the cutoff is removed. In other words, we 
demand that 

lim -^-GW(A, MA)) = 0. (23) 

The function /u(A) thus plays the role of the renormalized coupling constant. 

The dependence of the coupling constant on the cutoff is usually expressed in terms of 
the beta function, which is defined by 

M - A* (24) 

Within our formalism, Eqs. (|23|) and ([24]) can be used to calculate the beta function. 

Note that once Eq. ([23]) is satisfied and /x(A) is determined, then H e g (via Z), based on 
Crfcfc'(A, fj,(A)), should also be independent of A as A — > oo. Thus, the complete problem of 
renormalization is solved. 

We shall now illustrate the method prescribed above in two simple cases of a Dirac 
particle in 1 dimension and a Schrodinger particle in 2 dimensions [|3],|S| . In both these cases 
the interaction potential will be taken as a delta function in position space : 

V(x) = -fi 5^(x), (25) 

where n is the dimension of configuration space. In the momentum space the interaction 
potential would simply be a constant, i.e., 

V(k) = -/i . (26) 

We will choose H to be the pure kinetic operator, and our model space to consist of all 
states with momenta less than A. Thus Q projects onto the momentum range [A, oo]. 

With the choice of the interaction potential described above, the series in Eq. fl2"2| ) can 
be summed exactly and is given by 



Gkk' = t— % -J (k-k% (27) 

where I(u) is given by 

r°° 1 

I(u) = / cTp z— -. (28) 

v ; Jx y u-E {p) v ; 

Following the preceeding discussion we now introduce an ultraviolet cutoff A. Replacing fi Q 
by the renormalized coupling constant fi and using Eqs. ( p3|) and (plD , we obtain the beta 
function as 

/W = ^ A qX- ( 29 ) 

To obtain the explicit expression for the beta function we need to evaluate the integral 
appearing in Eq. (^H). For the 1 dimensional Dirac particle we have n = 1, E (p) = p + m 
and 

J(w) = / rfp ? r = - In ) x ( . 30 

The corresponding beta function is given by 

0(AO = -M 2 - (31) 

For the Schrodinger particle in 2 dimensions we have n = 2 and Eo(p) = p 2 /2 (we set 
the mass of the particle as 1.) Proceeding exactly as before, we obtain 

I^) = -2n\n(^pj (32) 

and 

(3 = -47r/i 2 . (33) 

Note that the results in both examples above have the desirable property that the beta 
function is independent of the model space cutoff, A. The beta functions calculated give rise 
to asymptotically free theories and generate the generally accepted pattern for the flow of 
the coupling constant for the two examples described above. 
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In conclusion, we have introduced a way of implementing a nonperturbative renormaliza- 
tion scheme within the context of many-body effective Hamiltonian theory. We have tested 
the method with applications to simple scale-invariant quantum mechanical systems. These 
examples exhibit the efficacy of our ideas and calculations for more realistic many-body 
systems and for quantum field theory are presently under investigation. 
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